Goldstone-type fluctuations and their implications for the amorphous solid state 
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In sufficiently high spatial dimensions, the formation of the amorphous (i.e. random) solid state of 
matter, e.g., upon sufficent crosslinking of a macromolecular fluid, involves particle localization and, 
concommitantly, the spontaneous breakdown of the (global, continuous) symmetry of translations. 
Correspondingly, the state supports Goldstone-type low energy, long wave-length fluctuations, the 
structure and implications of which are identified and explored from the perspective of an appro- 
priate replica field theory. In terms of this replica perspective, the lost symmetry is that of relative 
translations of the replicas; common translations remain as intact symmetries, reflecting the statis- 
tical homogeneity of the amorphous solid state. What emerges is a picture of the Goldstone-type 
fluctuations of the amorphous solid state as shear deformations of an elastic medium, along with a 
derivation of the shear modulus and the elastic free energy of the state. The consequences of these 
fluctuations — which dominate deep inside the amorphous solid state — for the order parameter of 
the amorphous solid state are ascertained and interpreted in terms of their impact on the statistical 
distribution of localization lengths, a central diagnostic of the the state. The correlations of these 
order parameter fluctuations are also determined, and are shown to contain information concerning 
further diagnostics of the amorphous solid state, such as spatial correlations in the statistics of the 
localization characteristics. Special attention is paid to the properties of the amorphous solid state 
in two spatial dimensions, for which it is shown that Goldstone-type fluctuations destroy particle 
localization, the order parameter is driven to zero, and power-law order-parameter correlations hold. 

PACS numbers: 61.43.-j, 82.70.Gg, 64.60.Ak 



I. INTRODUCTION 

The aim of this Paper is to identify the long 
wave-length, low energy fluctuations of the amorphous 
(i.e. random) solid state, and to investigate their physi- 
cal consequences. In particular, by constructing an effec- 
tive free energy that governs these Goldstone-type fluc- 
tuations, we shall determine the elastic properties of the 
amorphous solid, including its static shear modulus. We 
shall also analyze the effect of these fluctuations on the 
amorphous solid order parameter, and hence determine 
their impact on physical quantities such as the distri- 
bution of localization lengths and the order parameter 
correlations. Along the way we shall reveal the physical 
information encoded in these correlations. The treatment 
presented here is valid deep inside the amorphous solid 
state, i.e., far from the critical point associated with the 
phase transition from the liquid to this state. Fluctuation 
phenomena such as those discussed in the present Paper 
cannot be captured by a solely percolative approach to 
amorphous solidification. A brief account of the work re- 
ported in the present Paper has been given elsewhere . 

We shall pay particular attention to systems of spa- 
tial dimension two, for which we shall see that the ef- 
fect of fluctuations is strong: particle localization is de- 
stroyed, the order parameter is driven to zero, and order- 
parameter correlations decay as a power law in the sep- 
aration between points in the sample. Thus we shall 
see that the amorphous solid state is, in many respects, 
similar to other states of matter exhibiting (or nearly 



exhibiting) spontaneously broken continuous symmetry. 

This Paper is organized as follows. In Sec.[n]we sketch 
the properties of the amorphous solid order parameter, 
including the form it takes in the amorphous solid state, 
focusing on symmetry properties and how they mani- 
fest themselves within the replica formalism. In Sec. IIIII 
we describe the structure of the low energy, long wave- 
length Goldstone-type excitations of the amorphous solid 
state in terms of distortions of the value of the order pa- 
rameter. Here, we also make the identification of these 
order-parameter distortions as local displacements of the 
amorphous solid. In Sees. HVl and IVl we determine the 
energetics of these Goldstone-type excitations by begin- 
ning with a Landau-type free energy expressed in terms 
of the amorphous solid order parameter and ending with 
elasticity theory. Along the way, we derive a formula for 
the elastic shear modulus, which shows how this mod- 
ulus vanishes, as the liquid state is approached, at the 
classical (i.e. mean-field theory) level. In Sec. I VII we dis- 
cuss the impact of Goldstone-type fluctuations on the 
structure of the amorphous solid state by examining how 
they diminish the order parameter and modify the dis- 
tribution of localization lengths. We analyze the impact 
of such fluctuations on the order-parameter correlations, 
and also catalog the various length-scales that feature in 
the Paper. In Sec. lVIII we take a closer look at the effects 
of Goldstone-type fluctuations on structure and corre- 
lations in two-dimensional amorphous solids. In partic- 
ular, we show that Goldstone-type fluctuations destroy 
particle localization, the order parameter is driven to 
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zero, and power-law order-parameter correlations hold, 
and we illustrate these features for certain special cases. 
In Sec. IVIIII we discuss the physical content of order- 
parameter correlations in terms of spatial correlations 
in the statistics of the localization lengths. We also in- 
troduce distributions of correlators, and relate their mo- 
ments to order-parameter correlators. Some concluding 
remarks are given in Sec. IIXI Technical details are rele- 
gated to four appendices. 



II. AMORPHOUS SOLID ORDER 
PARAMETER, SYMMETRIES AND SYMMETRY 
BREAKING 



A. Order parameter 

The solid state of randomly and permanently con- 
strained matter, well exemplified by the rubbery state of 
vulcanized macromolecular matter, may be detected and 
characterized via the following order parameter, which 
depends on an arbitrary number v (> 2) of tunable (but 
nonzero) wave vectors {k 1 , k 2 , . . . , k^}: 



J 

3 = 1 



(i) 

Here, the vectors {Rj}j—i give the positions of the J 
particles that constitute the system; they inhabit a large, 
Z)-dimensional, hyper-cubic region V of volume V. An- 
gular brackets (• • • ) indicate an equilibrium expectation 
value taken in the presence of a given realization of 
the quenched random constraints, possibly in a state 
with spontaneously broken symmetry. Square brackets 
[• • • ] denote an average over the number and specifi- 
cations of the quenched random constraints. Periodic 
boundary conditions discretize the wave vectors k to the 
values 2irm/V 1 / D in terms of D-tuples of integers m. 
Additional discussion of this circle of ideas is given in 
Refs. 0,0. 

The order parameter detects and diagnoses the 
amorphous solid state by sensing and quantifying the 
presence of static random waves in the particle density, 
as implied by (e ik Rj ) ^ (for k ^ 0). To acquire some 
feeling for how it works, consider the illustrative example 
in which a fraction 1 — Q of the particles are unreal- 
ized whilst the remaining fraction Q are localized, har- 
monically and isotropically but randomly, having random 
mean positions (Rj) and random mean-square displace- 
ments from those positions 

((R, - <R,)) d (R, - (R,)),,) = S dd , (2) 

where d and d! are cartesian indices running from 1 to 
D. It is straightforward to see that for this example the 



order parameter becomes 
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is the disorder-averaged distribution of squared local- 
ization lengths £ 2 of the localized fraction of parti- 
cles HHEf 

The illustrative example correctly captures the pattern 
in which symmetry is spontaneously broken when there 
are enough random constraints to produce the amor- 
phous solid state: microscopically, random localization 
fully eliminates translational symmetry; but macroscop- 
ically this elimination is not evident. Owing to the ab- 
sence of any residual symmetry, such as the discrete 
translational symmetry of crystallinity, all macroscopic 
observables are those of a translationally invariant sys- 
tem. This shows up as the vanishing of the order param- 
eter, even in the amorphous solid state, unless the wave 
vectors sum to zero, i.e., Y2a=i k a = 0- 

The case v — 1 is excluded from the list of order param- 
eter components shown in Eq. Q). This case corresponds 
to macroscopic density fluctuations, and these are as- 
sumed to remain small and stable (i.e. non-critical) near 
the amorphous solidification transition, being suppressed 
by forces, such as the excluded-volume interaction, that 
tend to maintain homogeneity. Additional insight into 
the nature of the constraint-induced instability of the 
liquid state and its resolution (in terms of the formation 
of the amorphous solid state — a mechanism for evading 
macroscopic density fluctuations) is given in Ref. |3( , es- 
pecially Sec. 4.2. 



B. Symmetries and symmetry breaking; replica 
formulation 

How does the order parameter transform under 
translations of the particles? If, in the clement (exp ik a ■ 
Rj), one makes the translation Rj — > Rj + r a then the 
element is multiplied by a factor exp ik a -r a and so the or- 
der parameter acquires a factor expi k a ' ra - Now, 
in the fluid state no particles are localized and the or- 
der parameter has the value zero, so it is invariant un- 
der the aforementioned translations. By contrast, in the 
amorphous solid state the order parameter is nonzero, 
provided the wave vectors sum to zero. Thus, the or- 
der parameter varies under the translations unless the 
translation is common to each element, i.e., r a = r. To 
summarize, the liquid-state symmetry of the independent 
translations of the elements is broken down, at the amor- 
phous solidification transition, to the residual symmetry 
of the common translation of the elements. 



When replicas are employed to perform the average 
over the quenched disorder (i.e. the number and location 
of the constraints), what emerges is a theory of a field 
O(x) defined over (1 + n)-fold replicated space x, so that 
the argument x means the collection of (1 + n) position 
D-vectors {x°, x 1 , . . . , x™} conjugate to the wave vectors 
{k^k 1 , . . . , k"}. It is understood that the limit n — > 
is to be taken at the end of any calculation. In terms 
of replicas, the expectation value of this field (Q(x)) is 
proportional to 

and its Fourier transform fi(fc) = Jdx exp(ifc • x) Q{x), 
has the form 

The field Vl(x) fluctuates subject to the demand, men- 
tioned above, that the critical freedoms are only the 
corresponding Fourier amplitudes f2(fc) for which at 
least two of the D-component entries in the argu- 
ment k = {k°, k 1 , . . . , k"} are nonzero. Both a semi- 
microscopic approach and arguments based on sym- 
metries and length-scales yield a Landau- Wilson effec- 
tive Hamiltonian governing the fluctuations of this field, 
which is invariant under independent translations of the 
replicas but whose precise structure we shall discuss later. 
For now, let us just mention that the corresponding ex- 
pectation value of this field is the order parameter 
with v = 1 + n: it becomes nonzero in the amorphous 
solid state and, in doing so, realizes the pattern of spon- 
taneous symmetry breaking described above. Invariance 
under independent translations of the replicas breaks 
down to invariance under the subgroup of common trans- 
lations of the replicas. 

III. GOLDSTONE FLUCTUATIONS: 
STRUCTURE AND IDENTIFICATION 

A. Formal construction of Goldstone fluctuations 

In the amorphous solid state, one of the symmetry- 
related family of classical values of the order parameter 
has the form 

^e^WW, (6) 

in which W is real and depends only on the magnitude 
of k T . Some geometry is needed to define the variables in 
this formula. We introduce a complete orthonormal basis 
set in replica space {e a }" =0 i m terms of which vectors k 
are expressed as 

fc = V" k a e a . (7) 
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FIG. 1: A classical state in replicated real space: a hill in x 
space with its ridge aligned in the x\ direction and passing 
through the origin. The thickness of the lines is intended 
to suggest the amplitude of the order parameter, the thicker 
the line the larger the amplitude. Symmetry-related classical 
states follow from rigid displacements of the hill perpendicular 
to the ridge, i.e., in the x T direction. 

We also introduce the replica body- diagonal unit vector 

e=^_V e a , (8) 

relative to which we may decompose vectors k into lon- 
gitudinal (A) and transverse (r) components: 

k = k\ + k T , k\ = (k ■ e) e, k T = k — (k ■ e) e. (9) 

We find it convenient to parametrize the longitudinal 
components of position and wave vectors in the following 
distinct ways: 

1 - 

x x = (l + n) 1/2 x cm e, x cra = — — Vx Q ,(10a) 

1 + n ■* — ' 

a=0 

n 

kx = (l + nj-^ktate, k tot E^k a . (10b) 

Then x cm and k to t are, respectively, the analogs of the 
following conjugate pair of vectors: the center-of-mass 
position and the total momentum. With them one then 
has, e.g., k x - x x = k tot • x cm . 

The variables just introduced exhibit the structure of a 
classical state Eq. © as a rectilinear hill in a;-space, with 
contours of constant height oriented along xx, as shown in 
Fig. 2] The peak height of the ridge determines the frac- 
tion of localized particles; the decay of the height in the 
direction x T determines the distribution of localization 
lengths. The width of the hill corresponds to the typical 
value of the localization length. Symmetry-related clas- 
sical states are generated from Eq. by translating the 
hill rigidly, perpendicular to the ridge-line (i.e. parallel to 
x T ). Such a transformation corresponds to relative (but 
not common) translations of the replicas. 
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FIG. 2: Goldstone-distorted state in replicated real space: 
the hill is displaced perpendicular to the ridge to an extent 
Ut that varies with position x\ along the ridge. Note that the 
scale of this figure is much larger than that used for Fig.Q the 
thick line lies along the ridge of the classical state, but now 
it is the width of the line that indicates the width of the hill 
£t yP . The Goldstone-type fluctuations occur on wave-lengths 
longer than £ tyP - 

This pattern of symmetry breaking suggests that the 
Goldstone excitations of a classical state are constructed 
from it via x\ (or equivalently x cm ) -dependent transla- 
tions of the hill in the x T direction, i.e., ripples of the hill 
and its ridge. In two equivalent realizations, this gives 

VQ(k) = [ dx CB1 e ik ^- x '^ +ik ^ u ^ x '^W{k T ) 7 {llsi) 
Jv 

Va{x)=W(x T -u T {x cin )). (lib) 

Note that we are choosing to define Fourier transforms 
as follows: 

A(x) = J dk T tfktot e~ iK ^ e ^ ktot - x -" A(k), (12a) 
A{k) = Jdx T dx cm e lK ^e lkt ° t X ™l(x); (12b) 
also note that 

W(x T ) = Jdk T e~ lk T^W(k T ). (12c) 

Here and elsewhere, bars indicate division by factors 
of 27r; on integration measures there is one such di- 
vision for each variable of integration. The details of 
the excitation are encoded in the replica-transverse field 
M r (x cm ), an nZ?-component field that depends on the 
replica- longitudinal position x cm ; these are the Goldstone 
bosons, or phonon excitations, of the amorphous solid 
state. Provided the Fourier content of the Goldstone field 
M T (x cm ) occurs at wave-lengths long compared with the 
hill width (i.e. the typical localization length) — a point 
that we shall amplify below — these exhaust the spectrum 
of low energy excitations. 



The Goldstone excitations that we have just con- 
structed are analogs of the capillary excitations of the 
interface between coexisting liquid and gas states; see 
Ref. Q for a review. In the liquid-gas context they simi- 
larly accompany a spontaneous breaking of translational 
symmetry associated with the choice of interface location. 

What about other excitations of the broken-symmetry 
state? The complete spectrum of exciations is accounted 
for by decorating the Goldstone-type parametriza- 
tion (|llafl with additional freedoms w(k r ), so that the 
field £l(k) is expressed as 

VO(fc) = J dx cm e lkt ° t - x "» +lfcT - UT(x ™ ) (>V(A: T ) +w(fc)), 

. (13) 

where w is a real- valued field that depends on fc T , suitably 
constrained to be independent of the Goldstone excita- 
tions. 

The structure of the Goldstone excitations is in ac- 
cordance with the linear stability analysis of the clas- 
sical broken-symmetry state, due to Castillo et al. [||, 
which identified a family of linearly additive Goldstone- 
type normal modes of excitation indexed by k to t : 

VSl(k) = [ dx cm e lktot - x ™ W(k r ) + ik T -v r (k tot )W{k T ), 
Jv 

(14) 

with arbitrary replica-transverse amplitude u r (k to t)- 
That this linear glimpse of the Goldstone-type excita- 
tions is in accord with the nonlinear view focused on in 
the present Paper follows by expanding Eq. I|13|) to linear 
order in the Goldstone fields u T and omitting the non- 
Goldstone fields w, thus arriving at 

Vn(Jfe)M / dx cm e iktot X "»>V(fc r ) 
Jv 

+ik T - / dx cm e lktot - x ™u T (x cm )W(A: T ), (15) 
Jv 

which shows that u T (k to t) is the Fourier transform of 
M T (x cm ). The reality of f2(x) ensures that f2(— k) = 
fi(fe)* and, via Eq. ifltjl. that v T (— k to t) = v T (k tot )*. Via 
Eq. i|15|) . this in turn ensures the reality of u T (x cm ). 

Some insight into the structure of the Goldstone exci- 
tations, which induce deformations of the classical state, 
is obtained from its replicated real-space version. Con- 
sider the interpretation of f2(x) as a quantity propor- 
tional to the probability density for the positions of the 
1 + n replicas of a particle to have the values {x a }™ =0 ; see 
Eq. Then classical states, Eq. (|llb|> at constant it T , 
are ones that describe a translationally invariant bound 
states: fl(x) does not depend on the mean location of 
the replicas x cm , but does depend on their relative loca- 
tions, through x T , and decays the more the replicas are 
separated. This point is exemplified by the particular 
form given in Eq. I|26al) . Now, the Goldstone-distorted 
state, Eq. Ijllbfl with u T varying with x cm , also describes 
bound states of the replicas, but ones in which the de- 
pendence of the probability density on relative locations 
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FIG. 3: Molecular bound state view of the classical and 
Goldstone-distorted states in replicated real space. The 
full circles within a border represent the replicas of a given 
monomer location. Repetitions of them along a row indicate 
that the center of mass of the bound states is distributed ho- 
mogeneously. Upper bars: one classical state; the probability 
density is peaked at the shown configurations, in a manner in- 
dependent of the location of the center of mass of the bound 
state of the replicated particles. Middle bars: another classi- 
cal state, obtained by the former one via a relative translation 
of the replicas that does not vary with the location of the cen- 
ter of mass. Lower bars: a Goldstone-distorted state; the the 
probability density is peaked in a manner that varies with the 
location of the center of mass. 

varies with x cm . The particular form 126al) . which gives 

/>oo 

Jo 

xexp(-|.T r -M T (x cm )| 2 /2^ 2 ), (16) 

exemplifies this point; in particular, one sees that the 
most probable value u T of the relative locations x T now 
depends on the center- of-mass location x cm . These re- 
marks are amplified in Fig. [3J 

Returning to the issue of the structure and proper- 
ties of the Goldstone-type excitations of the amorphous 
solid state, recall that the critical Fourier amplitudes of 
the field are those that reside in the higher-replica sector 
[i.e. HRS, for which at least two D- vector elements of the 
argument of Sl(k) are nonzero]. Do the proposed Gold- 
stone distortions of the classical state excite the lower- 
replica sectors? If so, they would be suppressed by inter- 
actions, such as particle repulsion, that tend to preserve 
homogeneity. To see that they do not, let us examine 
the Goldstone-distorted state in the zero- (i.e. k = 0) 
and one- (i.e. k = qe Q with q ^ 0) replica sectors. In 
the former, one readily sees from Eq. (|lla|) that it has its 
undistorted value. A straightforward calculation shows 
that in the latter sector the distorted order parameter 
contains the factor: 

/ — jt- exp(iq-x cm + iq-u Q (x cm )), (17) 



where u° = u T ■ e a . By introducing the transformation 
x cm — * x' cm = x cm + u Q (x cm ) we see that the Goldstone- 
distorted state remains zero in the one-replica sector, 
provided all D- vector elements u" of the Goldstone field 
w T (x cm ) obey the condition 

\det(6 dd ,+d d u%)\ = l, (18) 

which, for small amplitude distortions, reduces to 
d d u d = 0. That this condition corresponds to in- 
compressibility (i.e. pure shear) will be established in 
Sec. 1111 Bl (Here and elsewhere, summations from 1 to 
D are implied over repeated Cartesian indices, such as 
d and d 1 .) This is as it should be: to make density fluc- 
tuations there must be some compression amongst the 
elements u a of u T . 



B. Identifying the Goldstone fluctuations as local 
displacements 

In this subsection our aim is to establish the connec- 
tion between the Goldstone fields and the displacement 
fields of conventional elasticity theory 0,0- Inter alia, 
this identifies the stiffness associated with the Goldstone 
fluctuations as the elastic modulus governing shear de- 
formations of the amorphous solid. As we shall see, the 
discussion is general enough to apply to any amorphous 
solid that breaks translational symmetry microscopically 
but preserves it macroscopically, in regimes where the 
constituent particles are strongly localized in position. 

Recall the amorphous solid order parameter, Eq. QJ, 

J- 1 V 3 (e lkl - R ^)(e lk2 - Rj ) • ■ ■ (e lk "- Rj )l , 
^-^j=l J 

and focus on a single element (e^'^). It is convenient 
to consider the element in the form 

e^p 3 (k), (19) 

where r, is the mean of the position of particle j, and 
p describes fluctuations about the mean. Now consider 
the impact of a long wave-length shear displacement, 
encoded in the field v(x), which deforms the probabil- 
ity density associated with the position Rj of particle 
j. By assuming that the typical localization length is 
much smaller than the length-scale associated with vari- 
ations of the deformation we are able to retain only the 
rigid displacement of the probability density for R^- and 
neglect any deformation of its shape. Thus, under the 
displacement the element is deformed as 

( e ik.R^ = e *k.r 3 p ^ k ) > e ik.(r 3 .+v(r 3 .)) ^ 

Inserting such deformations into the order parameter, 
with an independent displacement field for each element, 
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gives the distorted form 
1 J 

3=1 



p^k^-.-p^k") , (21) 



where we have arrived at the second form by noting 
that, in the amorphous solid state, the mean position 
of any particle is distributed homogeneously. Next, we 
decompose the collection of displacement fields {v°} into 
longitudinal and transverse parts, according to the geo- 
metrical prescription given in Sec. IIII Al but keeping in 
mind the fact that there are now v copies (rather than 
1 + n). The essential point is that the longitudinal part 
of {v a }, which corresponds to common deformations of 
the elements, does not generate a new value of the order 
parameter, in contrast with the transverse part (which 
generates relative deformations). Therefore, the physi- 
cal displacements are the transverse part of {v Q }. To 
see this point, consider the special situation in which the 
transverse part of {v a } is position dependent but the lon- 
gitudinal part is not. In this case, (V times) the order 
parameter becomes 

J dre lS «=i ka - r e^- t ' T(r)+lktot - v ™ [p^k 1 ) • • • p^k") 

. (22) 

As things stand, in the presence of v T (r), a longitudi- 
nal part v cm would have the effect of producing a new 
value of the order parameter. As, on physical grounds, 
we expect that it should not, we see that physical dis- 
placements are purely transverse. This is a consequence 
of the fact that, for the long wave-length displacements 
under consideration, the deformed amorphous solid state 
continues to preserve translational invariance macroscop- 
ically. This argument for the absence of the longitudinal 
part of the displacement field continues to hold when it 
has position-dependence: if, when constant, it does not 
generate a new state degenerate with the old one then, 
when varying, it should not generate a low-energy defor- 
mation. 

Thus we have realized the goal of this section: by com- 
paring Eqs. IjllaJI and ill' 21) we see that the formally- 
constructed Goldstone fields u T are in fact the physi- 
cal displacement fields v T . Actually, there is one fur- 
ther point to address, concerning the argument of the 
final factor in Eq. viz., [ p^k 1 ) • • • pj(k") ] . Ap- 

parently, there is dependence on the full set of wave 
vectors {k 1 , . . . , k"}, whereas Eq. (|llafl indicates depen- 
dence only on the transverse part of this collection. The 
resolution of this apparent discrepancy lies in the obser- 
vation that, for position- independent v T , Eq. (|22|l must 
revert to a classical state, for which the final factor does 
not depend on k tot . As the final factor consists of prob- 
ability clouds, which we are assuming to be undeformed 
by the displacements, this factor continues to be inde- 
pendent of k to t in Eq. (|22|l . 



IV. ENERGETICS OF GOLDSTONE 
FLUCTUATIONS; ELASTIC FREE ENERGY 

The simplest Landau- Wilson effective hamiltonian 
controlling the order parameter has the form 



S n = Vc J2 (-ar + ^ok-k)Q(k)fl(-k) (23) 

fcGHRS 
fcl,fe2,fc3£HRS 

where c is the number of entities being constrained per 
unit volume, and {ar, £o,g} are, respectively, the con- 
trol parameter for the density of constraints, the linear 
size of the underlying objects being linked, and the non- 
linear coupling constant controlling the strength with 
which the f2 fluctuations interact. HRS indicates that 
only wave vectors in the higher-replica sector are to be 
included in the summations. A semi- microscopic model 
of vulcanized macromolecular matter yields a = 1/2, 
t = (to 2 — m 2 )/m 2 , £ 2 = L£ p /2D and g = 1/6, where m 
controls the mean number of constraints (and has crit- 
ical value to c = 1), £ p is the persistence length of the 
macromolecules, and L/£ p is the number of segments per 
macromolecule . 

In terms of Sq, the replica partition function Z\ +n is 
given integrating over the critical modes of f2: 



Zi 



+ 71 



(24) 



Then, according to the replica scheme of Deam and Ed- 
wards the disorder-averaged free energy F is given by 



= lim n Zl+n v Zl = lim n JT^Zl+n, (25) 
1 n— >0 nZ\ n->o on 

where T is the temperature. 

It was shown in Ref. [3] (see also Refs. USDS) that 
making stationary with respect to f2 results in the 
classical state @), with 

/>oo 

W(fc T ) = Q / deme)e~ ek ' /2 , (26a) 
Jo 

Q = 2ar/3g, (26b) 



m 2 ) = {^l I are) ^1/ are), 



(26c) 



where tt{9) is the universal classical scaling function dis- 
cussed in Refs. 0,EJ. Distorting the classical state, via 
Eq. (|lla(l or Qllbjl . inserting the resulting state into Sq 
and computing the increase, S u , in iSn due to the dis- 
tortion u T (i.e. the elastic free energy) gives the following 
contribution, which arises solely from the quadratic "gra- 
dient" term in Sq: 

- TC ; 1+D I V n dk T ^W(k r ) 2 , (27b) 
(1 + n) 1+ ir J nD 



fin 
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in the former of which there are scalar products over both 
the nD independent components of u T and the D com- 
ponents of x. The derivation of this elastic free energy is 
given in App. [X] As we have discussed in Sec. IIII Bl and 
shall revisit in Sec. |VJ no is the elastic shear modulus. 
By using the specific classical form for W, Eq. I|26a|l . and 
passing to the replica limit, n — > 0, we obtain 



M0 - tcq 2 J dem 2 )de 2 m' 2 ) 

9 g 2 J e^ + e'- 1 



(28b) 



which, for the case of the semi-microscopic parameters 
stated shortly after Eq. <|23(l . becomes 



Ho = 2Tct 3 / dBdO 



(28c) 



The main technical steps of this derivation are given in 
App. [BJ Thus, we have arrived at the effective free en- 
ergy controlling elastic deformations, in the harmonic ap- 
proximation, ft is consistent with the result obtained in 
Refs. 11], in which the free energy cost of imposing a 
macroscopic shear deformation of the sample was deter- 
mined. 



V. IDENTIFICATION OF THE SHEAR 
MODULUS: MACROSCOPIC VIEW 



In Sec. IIII Bl we have explained why the Goldstone- 
type fluctuations are identified with local displacements 
of the amorphous solid by considering the impact of these 
fluctuations at a semi-microsopic level. In the present 
section we again address this identification, but now from 
a more macroscopic perspective, by coupling the Gold- 
stone fields to a force-density field. 

Accounting solely for the Goldstone-type fluctuations 
[i.e. ignoring the field w in the parametrization of the 
field O in Eq. (|13fl ], we approximate the replica partition 
function 1241) as 

Zi +n [f T ] ~ e~ Sn " l Jvu T exp ^ J dxd x ii T - d x u T 

+ 1^ dx/ T (x)-u T (x)), (29) 



in which Sn >c \ is the effective free energy (|23fl evaluated in 
the classical state ©, and T>u T indicates functional inte- 
gration over replicated displacement fields {u Q (x)} sub- 
ject to the following conditions: u is replica-transverse 
[cf. Eq. ©]; the D- vector elements it contains are pure 
shear [cf. Eq. 1|18|) ] ; and the Fourier content is restricted 
to wave-lengths longer than a short-distance cut-off £ < 
(which is we take to be on the order of the typical lo- 
calization length) but shorter than a long-distance cut- 
off t > (which is commonly on the order of the linear 
size of the sample). The reason for the restriction to 
wave-lengths longer than the typical localization length 
is that by restricting our attention to the Goldstone sec- 
tor of fluctuations we are omitting the effects of massive 
flutuations. To be consistent, we should also omit the 
effects of Goldstone-type fluctuations with wave-lengths 
sufficiently short that their energy scale is comparable 
to or larger than the scale for the (omitted) least mas- 
sive fluctuations. The appropriate criterion is that the 
short-distance cut-off be taken to be on the order of the 
typical localization length. This can be appreciated pic- 
torially from Figs. ^ and [3] Goldstone-type fluctuations 
having wave-lengths smaller that the hill width are omit- 
ted. Note that we have ignored the Jacobian factor con- 
nected with the change of functional integration variable 
from f2 to u. 

In order to reconfirm the identification of the shear 
modulus, we have, in Eq. I|29|l . coupled the displacement 
field u linearly to a replicated force density field f though 
a term — ^ J dx/(x) • u(x); because u\ is zero, only the 
transverse term, — i / dx/ r (x) ■ u r (x), remains. As for 
/ itself, it is taken to have D-vector elements that vanish 
in the zeroth replica and are identically equal to f in 
the remaining replicas. This reflects the fact that the 
force density is envisaged as being applied subsequent to 
the cross-linking process, and therefore does not feature 
in the replica that generates the cross-link distribution, 
but is repeated in the thermodynamic replicas (by which 
we mean the replicas that generate the logarithm of the 
partition function, not the disorder distribution). Thus, 
one has 



fr-fr = f-f-h-h = nt-t-(f-e) 2 



Tl re* n— >0 f. n 

f t — > ni l. 



1 + n 



(30) 



The integration over u in Eq. (|29|l is gaussian, and thus straightforward, requiring only the elastic correlator 
( u d{y) Ud'(y')), which is given in terms of the elastic Green function Gdd>{y — y') : 



, , , / pu exp (-J| J v dx d x u ■ <9 x u) u d (y) u d , (yQ T 

{ud(y)u d >{y )} = »— 7 - r — \ = — Gdd'{y-y ), 

J Vu exp (— J v «x a x u • a x uj ^ 



Gdd' (y) 



2-ir/ty 



'(Jke-^^k), 



G d d>(k) = (k 2 5 d d> - k d k d ,)/k A 



(31a) 
(31b) 



where £ > and are, respectively, the long- and short-distance cut-offs on the wave-vector integration, mentioned 
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above. This elastic Green function will be derived in App. [CJ In terms of the elastic Green function and the force 
density, one finds for the increase F[f] — F[0] in free energy due to the applied force-density field: 

F[f] - F[0] = -T lim Zl+ " [f] ~ r ^ 1+ " [0] = -T lim #- In (32a) 
LJ LJ n^o nZt[Q] n^odn Z 1+n [0] v ; 

« - — / dxdx' V / d (x)^(x-x')/ d '(x'), (32b) 
which indicates that /io is the shear modulus; see Refs. pll^- 



VI. EFFECT OF GOLDSTONE FLUCTUATIONS 
ON THE ORDER PARAMETER AND ITS 
CORRELATIONS 

In this section we discuss, how Goldstone fluctuations 
affect the expectation value and the correlations of the or- 
der parameter. Our discussion is based on harmonic elas- 
ticity theory as described by the free energy of ea. (|27a|) . 
We shall make frequent use of the elastic Green function, 
which is defined in eas. (|3rbfl and computed in App. [U] 
The effects of Goldstone fluctuations are most striking in 
two dimensions, hence the two-dimensional solid will be 
discussed separately in sec . l| Villi . 



A. Order parameter reduction due to Goldstone 
fluctuations 



Classically, the order parameter expectation value is 
given by Eq. 10 or, equivalcntly, 



(Vn(x,k T )) = ( /<2k tot e 



Vn(k)) (33e 



E 



(n(k)) « W(kr). (33b) 



The effect of Goldstone fluctuations on the order param- 
eter expectation value is estimated from the Gaussian 
theory (|27a|l . via which we compute the mean value of 
the distorted classical state: 



(vn(x,M) = ( /<?k to t 



-j'ktot-x 



VQ(k)) (34 



a 



e-**"*-* (fi(fc)) (34b) 

ktot 

ik r -u T (x) \ 



( 



] )yv(k T ). 



(34c) 



This is readily evaluated via the Gaussian property of u, 
and hence we find 

(Vn(x,k T )) w W(k T ), (35a) 
VV(fe T ) = exp(-TT D lg/2no) W{k T ), (35b) 

for the fluctuation-renormalized form of W(/c T ). Here 
D Tjj = ^<jd( x )|x=o and summation over repeated carte- 
sian indices is implied. Recalling the classical structure 



for W, Eq. I|26a|) . we see that the effect of the fluctuations 
is to induce, in say Eq. © or (|33b|) . the replacement 

poo 

W(k T ) = Q d£ 2 Af(t 2 ) e ' e ^ /2 — W(k T ), (36a) 



W(k T 



-Tr D kl/2fj, 



Q / deM{e-{TT D /^))e-^' 2 



JQ 



(36b) 



i.e. , a rigid shift of the distribution Af to longer (squared) 
localization lengths, as encoded in the new distribution 
Af. This is to be expected: the locally fluctuating local- 
ized objects are also subject to collective fluctuations — 
phonons. 

This shift of Af is determined by the value of which, 
as shown in App. El has the following leading-order de- 
pendence on D: 



L D 



D ~ 1 jn=>Z, for D>2; 

(37) 

for D = 2. 



(2ir) 2 D(D-2) Jl 



< • 



where S_d is the area of the (D — l)-dimensional surface 
of a D-dimensional sphere of radius unity, viz, 



Ez, = 2n D / 2 /T(D/2), 



(38) 



in which T(D), [not to be confused with T D of Eq. l|3*7|l ] 
is the conventional Gamma function. We see that in 
dimension D greater than two the shift is finite. How- 
ever, at and below D = 2 the shift diverges with the 
long-distance cut-off ^>, viz., the linear size of the sys- 
tem, doing so logarithmically in two dimensions. Not 
surprisingly, fluctuations destroy particle localization in 
two-dimensional amorphous solids and thus restore the 
symmetry broken at the classical level (see below). 

We now consider the scaling behavior of the 
fluctuation-induced shift of the distribution from Af to 
Af. We see from Eqs. (|36|l tha this shift is parametrized 
by the length £3, which is defined via 



£| = T d T/hq. 



(39) 
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How does this length compare with the typical localiza- 
tion length £typ; as defined, say, via the most proba- 
ble value of £ predicted by the classical theory? From 
Eq. I(26c() we see that 

£t 2 yP ~ e /r, (40) 

and thus we have that 

ejg yP ~r D TT/^e - (41) 

Now, from Eq. I(28c(l we have that 

Mo /T~ct 3 , (42) 

and we use this to eliminate fio/T in Eq. ijSJ. Further- 
more, from Eq. ((37(1 we have that, for D > 2, 

T D ~t%- D ~&- D , (43) 



provided we take the short-distance cut-off to be of order 
£ typ , as discussed in Sec.[V] We use this to eliminate T D 
in favor of £typ- Finally, by using Eq. 1)40(1 to eliminate r 
in favor of £typ/£o we obtain, for D > 2, 

^/e&p~(c& D )- 1 (w&) 6 " J '. ( 44 ) 

At D — 6 we find one-parameter scaling in the sense 
that £a ~ £typj up to the factor c£q , which measures the 
number of crosslinked entities within a region of order 
the size of a single one. This is to be expected, because 
we have used classical exponents for the divergence of 
£ty P and vanishing of /io with r, and six is the upper 
critical dimension for the transition to the amorphous 
solid state. In fact, if we were to replace these classical 
exponent by their anomalous values we would expect to 
recover one-parameter scaling at arbitrary D. 



I 

B. Two-field order parameter correlations 



As with their effect on the order parameter itself, the effect of Goldstone fluctuations on the two-field correlator 
can be determined from the Gaussian theory (|27a|l . which gives 

(VSl(x, k T ) FO(x', k' T )*) « ( e lfc ^(x) g -ifc> T (x')) w ^ (45) 

The required correlator (expi/c r • u T (x) exp — ik' T ■ u T (x')) is also readily evaluated via the gaussian property of u, 
and hence we have 



vn{x,k T )vn(x',k' T )*) «ex P 



\K - K\ 2 ^j exp (g dd ,(o) - g dd ,( x - x'))fc Td • k' Td ?j w(k T )w(k' T ), 

(46) 



2^0 

where the scalar products k rd ■ k' Td , and \k T — k' T \ 2 are respectively taken over n and nD components. 



r 



Recall, from Eq. (|37|l . that Te>\d = 2 diverges with the 
long-distance cut-off. This, together with the positive- 
semi-definiteness of the quadratic form (Gdd'(0) — 
Qdd'(^—^-'S)k T d-k' Td , in Eq. 1(46(1 . makes it evident that the 
correlator (V£l(x, k r ) VVt{-x! ', k' T )*) given in Eq. JJOJl van- 
ishes at D = 2 (and below) unless k — k'. This vanishing 
is a second facet of the fluctuation-induced restoration 
of symmetry discussed for the case of the order parame- 
ter, following Eq. [The correlator (Vfi(fc) VCl(k')*) 
vanishes in any dimension unless k\ — k' Xl owing to the 
preserved symmetry of common translations of the repli- 
cas, which encodes the homogeneity of the randomness in 
the amorphous solid state, i.e., its macroscopic transla- 
tional invariance.] If k = k! then, regardless of dimension, 
the correlator decays with increasing separation x — x', 
as shown by Eq. |@BJ. 

For D > 2 it is convenient to analyze the two-field 
correlator normalized by its disconnected part, i.e., 



Equations ((46(1 and 1(35(1 show this quotient to be given 
by 

expf— G dd .{r)k Td -k' Td , ) , (48) 
VMo J 

where the separation r = x — x'. 

We illustrate this behavior by considering the case of 
D = 3 and the regime £ < <C |r| <C £>, for which, as 
shown in Add. O (see also, e.g., Ref. Q), we may use 

G%(r)*-^(6 dd ' + r d r d ,), (49) 

where the unit vector f = r/|r|, and by choosing k = k' = 
{0, q, — q, 0, . . . , 0}. Then in this regime the normalized 
two-field correlator is given by 

exp (i^^r [1+cosV] )' (50) 

which depends strongly on the angle tp between the vec- 
tors q and r. 
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We may also examine the normalized two-field correla- 
te 

tor in the regime |r| <C £ < . Making use of Q dd i\p) m this 
regime, as given in Eq. icTTl of App.O we find that the 
normalized two-field correlator is given by 



/ |q| 

V37T/io £< 



/ 2T 



|q| 2 kl 2 



exp I — l t^ l ~t- [ — 2 + cos ip 



(51) 

where £ < = £ < /2tt (and t > = £ > /2ir). Of course, in 
this regime the result for the correlator is incomplete, as 
there will also be contributions from the non-Goldstonc 
excitations. 



C. Intermezzo on length-scales 

We pause to catalog the various length-scales featuring 
in the present Paper, and to indicate where they first 
appear. 

The shortest length is £ p , the persistence length of a 
macromolecule, which appears alongside L, the arclength 
of a macromolecule, shortly after Eq. (|23p. Together, 
they yield £o> the linear size of the objects being linked, 
which may be substantially larger than £ p . It is only 
through £o that £ p and L feature in the free energy, 
Eq. The density c of entities being constrained also 

first features in Eq. (|23|l . and sets a length-scale, c -1 / 15 , 
in Eq. I|44|) . Comparable to or longer than £o is the length 
£, the (distributed) localization length, first featuring in 
Eqs. J2Jl and ©. The typical value of £ is denoted £ t y P , 
first mentioned around Eq. I|40(l . The amount by which 
fluctuations shift the distribution of £ is encoded in the 
fluctuation length £a, first mentioned around Eq. (|39H . 
The short- and long-distance cut-offs for the Goldstone- 
type fluctuations are respectively denoted £ < and £ >1 and 
are first mentioned in Sec.[V] The linear size of the sam- 
ple is denoted £; see the beginning of Sec. IVIII 



VII. AMORPHOUS SOLIDS IN TWO 
DIMENSIONS 

We now focus on amorphous solids in two dimensions. 
By taking the long-distance cut-off £> to be the linear 
size of the sample £, so that its area is C 2 , we see from 
Eqs. H35fl and l|37(l that the expectation value of the order 
parameter does indeed vanish logarithmically with the 
sample size, doing with an exponent that depends, inter 
alia, on k T : 



(£ 2 ft(x,fc T )) « W(k T ) exp(-i77(^)ln(£/^<)) 



= W(k T ) (£/*<) 



-v(K)/2 



, (52) 



where the exponent r](n 2 ) varies continuously with wave- 
number and is defined via 



t?(k 2 



Tk 



(53) 



To arrive at this result we have made use of the two- 
dimensional real-space elastic Green function evaluated 
at the origin, computed in App.lCland stated in Eq. l(3"7)l. 
It confirms the expectation, mentioned above, that in 
two-dimensions fluctuations destroy particle localization 
and restore the broken symmetry. Note, however that, 
due to the purely entropic nature of the elasticity the 
exponent does not depend on temperature [cf. Eq. (|28c[) ]. 

Similarly, from Eq. I|46[) we see that the fc T -diagonal 
two-field correlator is given by 

(c 2 k T ) c 2 , k T y) « (54) 

W(k T ) 2 exp ( - — (G dd ,(0) - g dd >(r))k Td ■ k Td , 
V Mo 

where, as before, the separation is denoted by r = x — x' 
and the unit vector by f = r/|r|. In the regime £< <C 
|r| C, the Green function is given by [see Eq. (|C28|I in 
App.E] 

Gdd> (0) - Q dd > (r) = i- {6 dd , In (r/ £<) - f d f d , } , (55) 

where, for convenience, we have introduced the numeri- 
cally rescaled short-distance cut-off £ < , defined via 



(56) 



In consequence, the correlator decays algebraically with 
the magnitude r of the distance, and with an amplitude 
that depends on the orientation f of the distance (relative 
to k T ). Ignoring, for the moment the dependence on f, 
we have the leading-order result: 



(c 2 n(x, kr) c 2 n(x / , k T y) « w(K) 2 U</rj 



n(K) 



(57) 

As mentioned above, the fc r -off-diagonal two-field 
correlator — like the order parameter — vanishes in the 
thermodynamic limit, reflecting the fluctuation- induced 
restoration of symmetry. 

This scenario is similar to the theory of two- 
dimensional regular solids, taking into account harmonic 
phonons only, as discussed in Ref. The main differ- 
ence is the absence of Bragg peaks at reciprocal lattice 
vectors. Instead, there is a divergence in the scattering 
function [T^ at zero wave-number, as can be seen from 
the Fourier transform of the correlator (|57|) with respect 
to r (= x — x'): 

S(q,k T ) = J dre lr - q (£ 2 ft(x,fc r )£ 2 ft(x',fc r )*) 



W(k T ) 2 £ 



h(q£>) 




B(rj) («*<) 



»7-2 



B(x) = T 



xY{l - x/2) 
r(l + a;/2)' 



(58) 



In these formulae, Ji denotes a Bessel function and we 
have abbreviated n(k 2 ) as r\. One needs to keep the 
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two exhibited terms because there is an exchange of 
dominance as 77 passes through the value 1/2, i.e., at 
k 2 . = 2tt/io/T. Specifically, for (l>/l<) — > 00 at fixed 
q£ < the dimensionless factor {• • • } behaves as: 

{qi<)- 3/2 (1>/1<)?- V , for r/ < 1/2; (59a) 
(qt<) v ~ 2 , for 77 > 1/2. (59b) 

Restoring the dependence on orientation f , we arrive 
at the full, anisotropic form for the decay of the two-field 
correlator: 

£</r) exp f rj{kl) k Tdl ■ k Td2 f dl f d A ■ 

(60) 

Here, the symbol k Td indicates the unit vector k Td /\k T \. 
For the illustrative case of k = (0, q, — q, 0, . . .) this cor- 
relator reduces to 

/ r- \ 2 /~ \ 2r l{l 2 ) 

W{V2\ci\) (t < /r) exp (277^) cosV) , (61) 

where, once again, ip is the angle between q and r. 

The framework adopted in the present Paper allows 
the identification of the quasi-localized fraction Q, i.e., 
the fraction of particles that, whilst not truely localized, 
have R.M.S. displacements that diverge only logarithmi- 
cally, as the system size C goes to infinity. This fraction 
should be contrasted with the delocalised fraction, i.e., 
the fraction that have R.M.S. displacements that diverge 
linearly with C. One way to extract the quasi-localized 
fraction is via the order parameter l|52[) — specifically the 
amplitude of its power-law decay with C. A simpler route 
is via the two-field correlator (|6UII : 

Q 2 = lim lim (l < /r\ ^ (C 2 fi(x, k T ) L 2 0(x', k T Y) 

k r ^Q r—>oo \ J 

= W(0) 2 . (62) 

The option of analyzing the quasi-localized fraction 
would not be available to percolation-based approaches 
to vulcanized matter, as such approaches do not account 
for the thermal motion of the constituents (e.g. macro- 
molecules), but only the architecture of the structures 
they constitute. In fact, for the same reason, the entire 
circle of ideas described in the present Paper lie beyond 
the reach of percolation-based approaches. This shows up 
especially vividly in the context of the vulcanization tran- 
sition in low (and especially two) dimensions. Whereas 
percolation theory would indicate a lower critical dimen- 
sion of unity, with nothing fundamentally new happen- 
ing as one reduces the dimensionality through two, the 
present approach correctly finds a lower critical dimen- 
sion at two dimensions for the amorphous solidification 
transition. This is because, in addition to incorporating 
the physics of percolation [lj,[ll|, this approach contains 
the logically-independent physics of localization and the 
attendant issue of the spontaneous breaking of transla- 
tional symmetry. Thus, there are qualitative, and not 



merely quantitative, distinctions between the vulcaniza- 
tion transition and resulting quasi-amorphous solid state 
in two- and in higher-dimensional settings, owing to the 
strong role played by Goldstone-type fluctuations in re- 
duced dimensions. 

The results discussed in the present section echo 
those found for a variety of two-dimensional statistical- 
mechanical systems (for a survey see, e.g., Ref. 0), in 
the following sense. A sufficient density of random con- 
straints triggers a phase transition from a liquid state to a 
quasi-amorphous solid state. In the liquid state there are 
no clusters of constituents that span the system (i.e. the 
system does not percolate), there are no quasi- localized 
particles, order-parameter correlations decay exponen- 
tially with distance, and the static shear modulus is zero. 
In the solid state a cluster does span the system (i.e. per- 
colates), there are quasi-localized particles, the static 
shear modulus is non-zero, and order-parameter correla- 
tions decay alebraically, with a continuously-varying ex- 
ponent 77 that depends, inter alia, on the scale set by the 
probe wave-vector k T . However, the true long-range or- 
der found in higher dimensions, fails to set in, albeit only 
just, due to thermal fluctuations. 

The situation is reminiscent of that found in the set- 
ting of the melting of two-dimensional crystals, triggered 
not by the density of constraints, but rather by ther- 
mal fluctuations. In this setting, at high temperatures 
one has a non-rigid phase with exponentially decaying 
positional correlations associated with unbound disloca- 
tions. This is the analog of the non-percolating state of 
quasi-amorphous solids. At lower temperatures one has 
a rigid phase, associated with bound dislocations and al- 
gebraically decaying positional correlations, but no true 
long-range positional order. This is the analog of the per- 
colating regime of amorphous solids. The permanence of 
the architecture of the amorhphous solid forming sys- 
tems discussed in the present Paper prohibits the de- 
struction by thermal fluctuations of a percolating raft of 
constituents. However, percolation does not enforce true 
localization. 

In two-dimensional crystallization one also has the op- 
portunity for correlations in the orientations of the bonds 
connecting neighboring particles. Thus, at low temper- 
aures one has, in addition to quasi-long range positional 
correlations, true long-range orientational correlations. 
This brings the opportunity for an intermediate regime 
of temperatures, in which dislocations are unbound but 
disclinations, which would disrupt the orientational or- 
der, remain bound: the hexatic phase of two-dimensional 
crystals. 

This opportunity does not seem to be present in the 
model of vlucanized matter discussed here, which focuses 
exclusively on positional order and does not support 
topological excitations. Richer settings, such as those 
involving macromolecules with liquid crystalline degrees 
of freedom, seem likely to raise interesting opportunities 
for order in low-dimensional random systems. 
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VIII. PHYSICAL CONTENT OF CORRELATORS 



A. Identifying the statistical information in the two-field correlator 



What is the meaning of the correlators in the amorphous solid state of the vulcanization field theory? Let us begin 
with the two-field correlator (f2(^i) Q(k 2 )*}. Specializing to the case in which the zero-replica entries and k 2 are 
zero, the interpretation of this replica quantity is given by 



n(k 1 )Q(k 2 y 



J n 



7 £ n 



(63) 



jl,j2 = l "=1 



which expresses the connection with the (semi-microscopic) correlators of all pairs of particles that constitute the 
system. Consider a simple illustrative example, in which the positions of pairs of localized particles fluctuate about 
their mean positions according to a general gaussian correlated distribution. For such particles we would then have 

/gikj-Rjj -ik 2 -R J2 \ _ e ik ld r jld -ik 2d r j2de -±Ar jl j ldld2 k ldl k ld2 e -jArj 2 j 2Jldj k 2dl k 2d2 e Ar jlj2dld2 k ldl k 2d2 (64) 

where the mean, variances (J% — j 2 ) and co- variances (j% ^ j 2 ) of the positions are given by 



Tj = (Rj), 

= {{Rjxdi ~ (Rjidi)) {Rj 2 d 2 - (Rj 2 d 2 ))), 

and summations over repeated cartesian indices d\ and d 2 are implied. 

Putting such contributions together from all pairs of localized, particles we have 



(65a) 
(65b) 



= Q dri dr 2 dAru dAr 22 dAr i2 V(r 1, f 2 , Arn, Ar 22 , Ar i2 ) 



xe 



if idJ22=i k id~ if 2d'E2=i fc id e ~! A, "ii<ii<i2 Ea=i k id 1 k id 2e ~ 2 - Ar 22d 1 d 2 ES=i k 2d 1 k 2d 2 e ^r 12dld2 J22=i k Td 1 l 



(66) 



V(r 1, f 2 , Arn, Ar 22 , Ar i2 ) 



(QJ) 



-pp^2$(h- r n )6(r 2 - r j2 )5(Ar n ~ Ar njl )5(Ar 22 - Ar nn )8{Ar 12 - Ar hh ) 



ji,j 2 loc. 



where we have introduced V, the disorder-averaged joint distribution, over the pairs of localized particles, of the means, 
variances and co- variances of the particle positions. The macroscopic translational invariance of the amorphous solid 
state ensures that V depends on fi and f 2 only through their difference, and thus we replace V by V^^iii — 
f 2 , Arn, Ar 22 , Ar 12 ). By appealing to permutation symmetry including the zeroth replica, we can reinstate the 
dependence on the zeroth-replica wave vectors, and hence arrive at a hypothesized form for the two-point correlator: 

(fi(fci) n(*a)*) = Q 2 J dh dr 2 dAr n dAr 22 dAr 12 V' 1 V(h ~ f 2) Ar u , Ar 22 , Ar 12 ) 

Xe ifld E" = k ld~ir 2d ELo k ld e ~k Ar ^d 1 d 2 E" = k ld 1 k ld 2 e ~k Ar 22d 1 d 2 E^O k 2d 1 k 2d 2 g A »'l2d 1 d 2 ELo k Td 1 k 2d 2 _ (g7) 

What about contributions to the double sum over particles in Eq. (|63() associated with one or two unlocalized 
particles? Such contributions certainly exist, except in the limit of large crosslink densities (where all constituents 
are bound to the infinite cluster and are, therefore, localized); in the liquid state they would, of course, be the 
only contributions. There, they would give rise to diagonal contributions, i.e., the correlator (Q(ki) fi(fc 2 )*y would 
vanish unless k± = k 2 , owing to the intact symmetry of independent translations of the replicas. In the solid state, 
contributions associated with unlocalized particles arc expected to gives rise to short-ranged correlations. As such 
correlations are not the main focus of the present paper, we neglect contributions to the two-field correlator associated 
with unlocalized particles. 



B. Evaluating the statistical information in the two-field correlator 



In Sec. lVIIIAl we have connected the two-field correlator (V^(fc) Fil(fc')*) to the distribution V that characterizes 
pairs of localized particles [see Eq. 1)67(1] . and in Sec. I VI 51 have evaluated the Fourier transform of this correlator, 
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(Vil(x., k T ) Vft(x.', k' T )*), within the Goldstone-type fluctuation approach [see Eq. (|46ll ]. We now discuss the implica- 
tions of resulting correlator for the properties of the distribution. 

To do this, it is convenient to exchange the correlator in Eq. (|46|) for the Fourier transform: 



(il(k)Q(k')*) = 



dx. dx' 

V ~v~ 



e ik t0 fx e — ikJot- 



: ' (vn(x, fe r )vn(x',jfcj.)*) 



W(k T )W(K)5 ktot , K 



x exp ( — g dd > (x - x') k Td ■ k' Td , 



Mo 



T / rji 

— e tkt ° t x exp — ^(x) *Vd • fc; d / 



(68) 



where W(fc T ) is the fluctuation-renormalized order parameter given in Eqs. (|35bjl and l|36b[) . Next, in Eq. I|67[l we 
perform the integration over the center of mass of ri and r 2 , and equate the resulting form of the correlator to the 
form given in Eq. JfJHJ), having dropped the hats on the dummy variables, thus arriving at a formula obeyed by the 
distribution V: 



I gikitot-x 



Q 2 ^ t , k2tot /dxdA m dA. 22 dAn 2 ^x,Anr,Ar 22 ,An 2 )^ ^ 

xe "2 &r lldld2 'E2 = O k ld 1 k ld 2 g- 2 A ''22did 2 52a = o k 2d 1 k 2d 2 e ^ r 12dxd 2 Y,2 = Q k ld 1 k 2d 2 

= W(fc lT ) W{k 2T ) 5 kltot , k2tot J y e lk »- x exp (j- G dld2 (x) k lTdl ■ k 2Td .^j . (69) 

By solving this equation for T-^x, Arn, Ar 22 , Ari 2 ) one learns that it has weight only at values of (Arn, Ar 22 , Ari 2 ) 
of the form (£fl, (T/fi )Q dld2 (y)j for some values of the parameters (£i,£i,y), where 1 is the identity in D- 
dimensional cartesian space. Anticipating this, it is convenient to introduce the following parametrization in terms of 
a reduced distribution P(x, £f , £|, y): 

7>(x,Ar n ,Ar 22 ,Ar 12 ) = J dgN {£) d$Af(£) S(Ar u - 6(Ar 22 - $l) 

x J ^ * (Ar u - (T/ Mo ) 0(y)) P(x, y). (70) 

Note that ^(x, £f , £|, y) only becomes a true distribution when an appropriate Jacobian factor associated with the y 
dependence is introduced; nevertheless we shall continue to refer to it as a distribution. As shown in App.[Dl Eq. ltBU|) 
can be solved for the reduced probability distribution; the solution is given in Eq. (|D9() . Thus one finds for the full 
distribution: 

P(x,Ar n ,Ar 22 ,Ar 12 ) = J d&tf(g) dgfi fe 2 ) S(Ar n - gl) <5(Ar 22 - £ 2 2 l) 

X Jv{I ^ qe_lq(X ~ y) exp + q2 ~ Ar ^ ld2 Q dl fe) | S(Ar 12 - (T/^)Q(y)). (71) 



This has the structure of a source term, associated with the range of values of Q, convoluted with an appropriate 
"propagator" 

It should be borne in mind that, in this solution for the full distribution, Eq. |J7TJ, the integration over wave vectors 
q is subject to the usual cut-offs, i.e., those featuring in Eq. I|31b|) . Thus, even though the exponent in the factor 
exp (!(£? +£ 2 ) I 2 ~ Ari 2dld2 q dl q d2 ) grows at large q (i.e. the factor is not a decaying Gaussian factor), the large- 
wave- vector cut-off protects against divergence. In fact, if we focus on the distribution V at separations x that are 
large compared with the cut-off ^< then for values of £i and £ 2 with appreciable weight the aforementioned exponent 



14 



is small, and we may expand to obtain 

P( X ,Ar u ,Ar 22 ,Ari 2 ) = j d^Kr^l) d&N (£) 6(Ar n - S(Ar 22 - £ 2 2 l) 

expf-^ef+eDV^+Ar^^^^) J ^ £(x - y) <5(Ar 12 - (T/ Mo )0(y)) (72a) 

x ^ cx v(-l(£ + £)Vl + Ar 12dld2 d dl d d2 ^j $(Ar 12 - (T//*o)0(x)) (72b) 

« / ^ 2 2 -^ fe 2 ) «(Ar u - ^(Ar 22 - e 2 2 l) 

^ |l - i(ei 2 + CI) V 2 + Ar 12dld2 8 dl d d2 } S(Ar 12 - (T/ Mo )(?(x)) . (72c) 



x ■ 
V 



Here, d d = d/dx d , and S indicates the smoothed delta function resulting from the application of the wave-vector 
cut-off to the Fourier integration that yields it. We have, however, proceded to the second and third lines without 
explicitly indicating the effect of this smoothing, viz., the replacement of the factor #(Ari 2 — (T//Xo)^(x)) the the 
same quantity smeared over a region around the point x having linear dimension of order the short-distance cut-off. 



C. Two- (and higher-) field correlators as distributions of particle correlations 

Let us pause to revisit the interpretation of the two-field correlator. To do this, we introduce the disorder-averaged 
distribution Ai 2 of two-particle correlators C, viz., 



M 2 [c]=[r 2 ]T n <Mc klka -(e iki - R * +<k *- R «)) 

ji J2 = l ki,k 2 #0 



(73) 



in which 5 c (x + iy) = S(x) S(y). We can then observe that the two-field correlator, Eq. (|63fl . can be expressed as a 
suitable moment of M. 2 : 



/n 
VCM 2 [C] \{C KK 
^ — n 



(74) 



Note that we have appealed to replica permutation symmetry in order to reinstate the dependence on the zeroth-replica 
wave vectors. 

It is straightforward to extend this discussion to the general case of r-particle correlators 
^giki-Rjj+iki-RjjH hik r R jr ^ f or wn j cn the distribution M r is given by 



M r [C] 



r 1 
IF 



e n 



k 2 -k, 



_ / iki-R^+iki-Rj-jH Hk r -Rj 



(75) 



3i,...,jV=lki,k 2 ,...,k r ^0 

Observe that the r-field correlator can also be expressed as a suitable moment: 



« n 

<fi(Ai)n(Aa)-..fi(A P ))= VCM r [C] n c k f k°...k 



(76) 



Again we have appealed to replica permutation symmetry in order to reinstate the dependence on the zeroth-replica 
wave vectors. The distributions M. r are natural generalizations of the distribution of local density fluctuations, 
explored, e.g. in Ref. [Toj . 



IX. CONCLUDING REMARKS 

In this Paper we have identified the long wave-length, 
low energy Goldstone-type fluctuations of the amorphous 



solid state, and investigated their physical consequences. 
By constructing an effective free energy governing these 
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fluctuations, we have determined the elastic properties 
of the amorphous solid, including its static shear mod- 
ulus which, we have re-confirmed, vanishes as the third 
power of the amount by which the constraint density ex- 
ceeds its critical value (at the classical level). We have 
also analyzed the effect of these fluctuations on the amor- 
phous solid order parameter, finding that, in spatial di- 
mensions greater than two, they induce a simple, rigid 
shift of the distribution of (squared) localization lengths. 
In addition, we have explored the properties of the order- 
parameter correlations in the amorphous solid state, es- 
tablishing their physical content in terms of a joint prob- 
ability distribution characterizing pairs of localized par- 
ticles. Moreover, we have computed the corresponding 
correlator induced by Goldstone-type fluctuations and, 
hence, obtained a specific formula for this joint probabil- 
ity distribution. 

We have paid particular attention to systems of spa- 
tial dimension two. In this setting we have shown 
that fluctuations restore the symmetries broken spon- 
taneously at the classical level, particle localization is 
destroyed, the order parameter is driven to zero, and 
order-parameter correlations decay as a power-law in the 
separation between points in the sample. The state is 
a quasi-amorphous-solid state, inasmuch as it possesses 
algebraically-decaying correlations and rigidity. 

Our work can be extended in several directions. We 
have focussed here on the critical behavior of the shear 
modulus. It is also of interest to study elasticity in 
the strongly crosslinked limit, in which (geometrical and 
thermal) fluctuations are less important and mean-field 
theory should even be quantitatively correct. Whereas 
for long-chain macromolecules this limit is difficult to 
achieve, it may well be possible for Brownian particles, 
i.e., networks that are built from small units, in the sim- 
plest case just monomers [lfj. Even more interesting 
is a generalization to nonlinear elasticity Rubber can 
withstand very large deformations — up to 1000% — thus 
allowing large amounts of elastic energy to be stored in 
the system. Our approach can be generalized to include 
nonlinear terms in the strain tensor, which will arise from 
two sources: the expansion of the Goldstone type fluctua- 
tions (|llaf) and higher-order nonlinearities in the Landau- 
Wilson free energy (|23|1 . Work along these lines is in 
progress. 
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APPENDIX A: FROM LANDAU- WILSON 
HAMILTONIAN TO ELASTIC FREE ENERGY 

In this appendix we show how to get from the Landau- 
Wilson effective hamiltonian Sq, Eq. 1)23(1 . to the elastic 
free energy (|27b|) . Specifically, we compute the increase 
S u in Sq when the classical value fi c i, Eqs. © and J5BJ|, 
is replaced by Goldstone-disorted classical state (lllafl . 
parametrized by u T (x). There are three terms in (|23f) 
to be computed. Two are "potential" terms, which we 
shall see to have no dependence on u T (x); the third is a 
"gradient" term, and this is the origin of the dependence 
of S u on u T (x) . 

Before focusing on any individual terms, we note that 
we can express summations over higher-replica sector 
wave vectors as unrestricted summations, less lower- 
replica sector contributions; e.g., 



E 



feeHRS 



E -E -E 



(Al) 



Here, the two components of the lower-replica sector, viz., 
the one- and zero-replica sectors, are respectively denoted 
IRS and 0RS. 



J 



Applying this sector decomposition to the first term in Eq. (12.''>[) , and recognizing that there is no contribution from 
the one-replica sector and that the zero-replica sector contribution is simple, we have 



V |fi(fc)| 2 = (V-V -V ) \n(k)\ 2 = Y \n(k)\ 2 -o-Q 2 



(A2) 



Recall that we are concerned with the increase in Sq due to the Goldstone distortion. Evidently, of the contributions 
considered so far only the unrestricted sum has the possibility of being sensitive to the distortion. However, as we 
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shall now see, not even this contribution has such sensitivity: 



ot -x 2 g~ *t ■ «t (xi )+ife T -u T (x 2 ) 



W(fc T ) 5 



^EW = y J Oxter J ^^e-*--+*- 

<ffc r W(fc T ) 2 e -ifer-K(x 1 )-«.(x 2 )) y ^ e -ik tot .( Xl -x 2 ) 

""dfcr W(fc T ) 2 e -t T 'K("i)-^(x 2 ))( 1 + n )-f 5 ( Xl _ X2 ) = (i + n )-f / <ifc r W(fc r ) 2 , (A3) 



= V 
= V 



rfxi <ix 2 

F V 



independent of u T (x). Note that here and elsewhere in the present Apendix we shall anticipate the taking of the 
replica limit by omiting factors of V n . 

Turning to the third term in (|23|) . the nonlinearity, and handling the constraints on the summation with care, we 
are faced with the term 



iy-iki+ki+ks) 



V ^2 8k 1 +k 2 +k 3 ,o n(fei) f2(A?2) fi(A;3) = J dkidk-zdka j dy< 

fci,A; 2 ,fc 3 

X J d*l rfx 2 dx 3 e 'ki to fXi+ t k 2tot .x 2 +ik 3tot .X3^ lT . tlT (x 1 )+.t 2T -« T (x 2 )+ ! fc3 T .« T (x 3 ) W(Al T ) W(*2r) W(fc 3r ) 

= J dk lx dk 2X dk 3X J ^e^ A (felA+fc2A+fc3A) J dk lT dk 2T dk 3r J dy T e^- (k ^ +k ^ +k3 ^ 

X f dxi dx 2 dx 3 e ^i t „ t -x 1 + l k 2tot .x 2 + ik3 tot -X3 e ifc 1 .-« T (x 1 )+ l fe2 T -« T (x 2 )+ l fc3 T -u T (x 3 ) yy^r) W(fc 2T ) W(fc 3 r) 

= (1 + n) D J dkidk a £ka J ^ e -*y(ki+ka+k») J dk lT dk 2r dk 3T J dy T e - iy -< k ^ +k2 - 



r+fe 3T ) 



/ dXl dx 2 dx 3 e * 1 -X 1 +ik 3 .X 3 +ik3-X3 e ife lT .« T (x 1 )+ife2r-«r(x 2 )+ife 3T -«r(x3) W(Jb lT ) W(fc 2 r) W(fc 3T ) 

(l + n) D Jdki T dk 2T dk 3T S klT +k 2T +k 3T ,o J J dx 1 dn 2 dx 3 6(xi - y) <5(x 2 - y) £(x 3 - y) 



Xfi 



ifci T -u T (xi)+zfc 2T -u T (x 2 )+ifc 3T -u T (x 3 ) 



W(fci T )W(fc 2r )W(fc 3T ) 



= (1 + n) D J dk lT dk 2T dk 3T 4 lT+fe2T+fc3T ,o / ^ e '(^+^+^)^Ww(fc lT )w(fe r )W(^) 

= (1 + n) D y dk lT dk 2T dk 3T 4 lT +fc 2T +fc3 T ,o W(fcir) W(/c 2t ) W(fc 3r )- 
This is independent of u r (x). 



(A4) 



r 



Turning to the second term in (|23|l . the gradient term, derivatives: 
we have 



^ k ■ k I fi ( fc )! 2 = ^ dk T (k X -k X + k T - k T ) 



Y, k ■ k l r! ( fc )l 2 = v [ dk x dk T k x ■ k x 

7, ** 



xe 



rfXi dx 2 ik tot . Xl+ fk tot .x 2 
ifc T -« T (xi)+ife T '« T (x 2 ) yy^^^ 2 



dxi cfx2 -iktot-xi+iktofXa 

V V 



xe 



-ifc T -« T (xi)+ife T -u T (x 2 ) 



W(fc T ) 5 



(A5) 



(l + n) _£,/2 (l + n) _1 y dk T dk 



Of the two contributions airising from this term, from fc 2 



-ifc T -u T (xi) gife T -« T (x 2 ) yy^ ^ 2 



ik tn fX 2 \ 

(A6) 



and from fc 2 , the latter has no u T (x) dependence. This xe 
follows via the mechanism that we saw for the first term, 
viz., the development of a factor of <5(xi — x 2 ). As for the 

former contribution, to evaluate it we replace k x • k x by Next, we integrate by parts, once with respect to Xi and 
(1 + n) _1 k to t • k to t and generate this factor via suitable once with respect to x 2 , to transfer the derivatives to the 
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exponentional factors containing u T , arriving at 



which proceeds as follows: 



J dktot e -*—( X1 - x ^ j dk T W(k T ) 2 
(d- ixl e- iK - u ^A ■ (9 lX2 e^- u ^ X2 )) . (A7) 



Powers of k\ higher than two would, via the correspond- 
ing derivatives, produce higher gradients of well as 
nonlinearities involving lower-order derivatives. By per- 
forming the derivatives that we have here, as well as the 
ktot integration, which generates a factor <5(xi — X2) and 
allows us to integrate over, say, X2, we obtain 



(1 + n)- 1 - ^- 1 J dk T W{k T f 

x / dx (fc T • <9 x u r (x)) • (k T ■ <9 x u r (x)) , (A8) 



where the scalar products inside the parentheses are over 
n£)-component replica-transverse vectors whilst the out- 
side scalar product is over _D-component position vectors. 
The next step is to observe that, owing to the rotation- 
ally invariant form of W(fc T ), Eq. (|26a[) . the k T integra- 
tion includes an isotropic average of two components of 
k T , and is therefore proportional to the identity in nD- 
dimensional replica-transverse space. Thus, we arrive at 
the form 



(l + n)- 1 - / 2 



nD 



3k T W(k T ) 2 k T ■ k T 



rfx 

V 



(d x u T (x)<9 x u r (x)) , (A9) 



which involves the two types of scalar product mentioned 
beneath Eq. (| A8|) . Reinstating the factor of Vc from 
Eq. I|23|) , we identify the stiffness divided by the temper- 
ature) fi n /T and, hence, arrive at Eqs. l(77|) . 



APPENDIX B: EVALUATING THE STIFFNESS 
(A.K.A. SHEAR MODULUS OR RIGIDITY) 

In this appendix we display the main steps for obtain- 
ing the stiffness /io, Eqs. (|28l) . by evaluating the formula 
for /i„, given in Eq. I|27b|) and derived in App. in 
terms of the elements of the classical state, Eqs. (12611 . 
This involves the evaluation of the RHS of Eq. (|27b(l . 



dkr k 2 T W(k T ) 2 

= q 2 J dk T k 2 JleM(e) d£' 2 jv (£ >2 ) e -(« 2 +« ,2 K/ 2 
- q 2 Jle M{e) di 12 N{e) J & T k 2 T e-^ 2 +^>'/ 2 

= Q 2 / ^ 2 AA(e 2 )< 2 AA(e' 2 ) 
Jo 



a dk T e- Ak '' 2 



xd -A/2\ A=i 2 +i ,2 

DO 

dc 2 AA(e 2 x 2 AA(r 2 ) 

xd. A/2 \ A=e+e . 2 (27rA)- nD / 2 
= Q 2 j d£ 2 A/"(£ 2 ) d£' 2 A/"(£' 2 ) 

<nD{2^A)- nD ' 2 A- 1 \ A=e+il2 



n— »0 



poo 

idq 2 / deN{e)dex(e) {e+ey 1 m 

Jo 



APPENDIX C: ELASTIC GREEN FUNCTION 

The elastic Green function Qdd'{^ — x') featuring in 
Eq. I|31a|) arises via functional integration over the dis- 
placement field u(x). This integration comprises fields 
configurations that are (i) volume-preserving (at least to 
leading order in the gradient of the displacement field), 
and (ii) have Fourier content only from wave-lengths ly- 
ing between the short- and long-distance cut-offs and 
l > . As the weight in Eq. i|31a|) is gaussian with respect 
to u(x), we have that (T//xo) <?<jd'(x — x') is proportional 
to the inverse of the operator appearing sandwiched be- 
tween two displacement fields in the exponent of the 
gaussian, 



Y J d-K (<9 x u • <9 x u) 



(CI) 



provided this inverse is the one associated with the 
Hilbert space of vector-field configurations contributing 
to the functional integral. An application of the diver- 
gence theorem shows that, in addition to conditions (i) 
and (ii), the Green function obeys the equation 



-V 2 <?dd'(x) = 6 d d> <5(x), 



(C2) 



in which the delta function is to be interpreted as the 
identity in the appropriate Hilbert space, mentioned 
above. To determine Odd' ( x ) we express it in its Fourier 
representation: 



GM*) = Jdke- lk *g dd 



(k). 



(C3) 
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Next, we insert this representation into Eq. (|C2|I to ob- 
tain 

-V 2 /^ke- ikx ^(k) = Jdkk 2 e-^G dd/ (k) 



2tt/ 



<fk e - lkx (<W ~k- 2 k d k d >). (C4) 



where t < = £ < /2tt and £> = £ > /2tt. The final term, the 
integral representation of the appropriate delta function, 
accommodates restrictions (i) and (ii). Then, from the 
linear independence of the plane waves, we see that the 
Fourier integral representation of the Green function is 
given by Eq. I|C3(1 . with the amplitude given by 



(fc 2 5 dd , - k d k d ,)/k 4 , for f > 1 < k < 



otherwise; 



(C5) 

as given in Eq. I|31b|) . 

Before computing the real-space form of this D- 
dimensional Green function, we evaluate it at argument 
x = 0, as well as at small argument (|x| <§; ^<). At x = 
we have 

£drf'(x)| x =o = fdk (fc 2 S dd , - k d k d >) fc" 4 (C6a) 

J h.c.o. 



= < W° — 1 / tfkk 2 



Tn = 



D 

D-l E D 



D (2tt) 



(C6b) 
(C6c) 

dkk D - 3 . (C6d) 



Here and elsewhere, the subscript h.c.o. indicates that the 
integration is subject to the hard cut-off shown explicitly 
in Eq. JC6dJ). Evaluating the last integral for D > 2, 
and for D > 2 retaining only the dominant contribution, 
gives T D , as given in Eq. (|3TJl . 

Now generalizing to |x| small but nonzero, i.e., |x| <C 
i?<, we have, by expanding the exponential in Eq. (|C3|) . 

£w(x) « 6 dd ,T D - -[(D + l)5 dd , -2x d x d ,]\x\ 2 



(2ir) D D(D + 2) 



dk k 



D~l 



(C7) 



<W T D --[(D + 1)<W - 2x d x d ,\ |x| 2 



^ D 2 (D + 2)' 



(C8) 



where, again, we have retained only the dominant con- 
tribution. Specializing to D = 2 and D — 3 we find 



g ( d 2 d \ (x) » <w r 2 - -f- ^- [3 <w - 2^ , (C9) 



l_w 2 

64tt f 2 



0$ (x) « <W T 3 - -V J^- [2 5,,, - ^ 4*] . (CIO) 



457r£< t\ 



We now compute the real-space form of this D- 
dimensional Green function, valid for arbitrary |x|: 

Gm>(x) = Jdke-^g dd ,(k) (cii) 

= / ^ke- ik x (k 2 S dd >-k d k d ,)k- A (C12) 

Jh.c.o. 

= -(<WV 2 -d d d d ,)H(x), (C13) 
H{x) = f (?ke- ik x fc" 4 (C14) 

Jh.c.o. 

p ^ p-ik^x 

(C15) 



cffcfc' 



D-l l d D-l^ 



fc 4 



Specializing to the case of D — 3, and using spherical 
polar co-ordinates (fc, 8, tp), we have 

r*< rib r 27r 

H (3) (x) = (2^)- 3 / ^ dtp d9 sin0 e- lk ^ cose 



(C16) 



o ./o 



|x| r\ x \/ e < gin ^ 
271 " 2 V|x|/f> z3 



To control the potential divergence at small z we add and 
subtract the small- z behavior of sinz, thus obtaining 

= K(si 3 (|x|/f < )~Si 3 (|x|/£ > ) 



2tt 2 ' 



Sis (t) — dz 



sm z — z 



(C18) 
(C19) 



where Si 3 (z) is a generalized sine integral, which has 
asymptotic behavior 



Si 3 (t) 



-i/3! 



for i < 1, 



-(tt/4) - t" 1 , for i > 1. 



(C20) 



Using this behavior to approximate H^(x) in Eq. (|C18|I . 
and inserting the result into Eq. I|C13|I . noting that the 
final term vanishes under differentiation, we obtain for 

l< < |x| < t> 



1 



8tt|x 



(<Jdd' + idid') (C22) 



Now specializing to the case of D — 2, and using plane 
polar co-ordinates [k,ip], we have 



H {2) {x) = (2tt)- 



i fc 



— / ^e^ fc l x l cosv 



|2 /•*, 



2tt 



|x|A> 



(C23) 
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To control the potential divergence at small z we add Noting that S(t) has asymptotic behavior 01 
and subtract the small- z behavior of the integrand, thus 

obtaining { 4,t 2 , for t « 1, 

«£(*) ~ \ ( C27 ) 
H ( 2)(x) = ]x£ ^(^.^^ \ihx(t/2e— ), for**!, 

2tt J\x\/t> z Jo 

I |2 /-|x|/£< ^ z where 7 (= 0.5772 . . .) is the Euler-Mascheroni constant 

+Y~ / —( 1_ l z2 ) ( C24 ) and using this to approximate H^(x) in Eg. (|U25) . and 

J l x l/ £ > inserting the result into Eq. 1|C13J) . noting that the con- 

Ul2 (s(|x|/£ < )-S(|x|/£ > )) + " 



f > — f < stant term vanishes under differentiation, we obtain for 



— ^yW^-^W^jj-r—j— £ < «|x|«£ > the result 



-1J- In (*>/*<). (C25) 

rt j_ /•27T 



^ " " ^(x)w-^ : (<S (i<i /ln(e^ + 3|x|/2f > )-i (1 ^). (C28) 



r dz 

S(t) ee / — / ^ ( e -» ^-[i- iz 2 ]). (C26) 
Jo z Jo 



J 



APPENDIX D: FROM THE CORRELATOR TO THE DISTRIBUTION 

In this appendix we give the technical steps involved in going from the two-field correlator to the distribution V, 
as discussed in Sec. IVIII Bl Beginning with Eq. I|69|) . setting k2tot to be knot, and reorganizing, we obtain 

x y'dxe lktot - x 7'(x,Ar 11 ,Ar 2 2,Ar 12 ) = W{k lT )W(k 2r ) /ye^ew ^ (,)tlT * ,4H> . (Dl) 

Inserting the paramctrization (|7U|I for V, performing the resulting integrations over Arn, Ar 22 and Ari 2 on the LHS, 
exchanging factors of W for A/" via Eq. I|36b[) . cancelling factors of Q, and identifying the Fourier transform 

^(q,Ci 2 ,Ci,y) = y^ xe ^ x nx,£i,£!,y), (D2) 

we arrive at the folloinwg equation: 

/ A^fe 2 ) rfyP(k ltot ,^,e 2 2 ,y) e ^ 5?fe? -^ fci+ - 

= / <%lti(g) ^-^fe 2 ) e-'^*^-*^*^ /" dye 4kltotX e^ 0-1^) ^-fer^ (D3) 

In order to clarify the content of this equation, we re-write the wave-vector dependence on the LHS in terms of 
replica- longitudinal (more precisely, k T ) and replica-transverse (in fact k to t) components [see Eqs. <|9ll(Jb[l ]. noting 
that kitot — k 2to t and writing q for each, to obtain 

Next, we identify Laplace transformations with respect to £f and £§, equate the entities being transformed on the 
LHS and RHS, and take the replica limit, thus arriving at 

e -H£ + £) c ' 2 J dyV(q, £i,£f,y) e« G dld2 (y)k lTdl -k 2rd2 e ^ Q dld2 {y) q dl q d2 _ J dy e itl y e^o Sd 1 d 2 (y)k lTil -k 2T d 2 _ ^g-j 

The next steps are to move the gaussian prefactor to the RHS, and to introduce the dummy variable g, which takes 
on the values held by the second-rank tensor Q (y) : 

= J dge^ **i*3 k ™i- h *r<h> e h(&+&y J dye lCiy 6(g-g{y)). (D6) 
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Again equating entites being transformed, this time being transformed with respect to g, we find 

J dy5(g-g(y))V(q,ei,e2,y) = eKS+SV f dye i ^J(g - G(y)) e"* ^^(y)^^. (D7) 

Next, we equate entites being transformed as J" dy 6 (g — £/(y)) ■ • • , thus arriving at the Fourier transform of the 
reduced distribution: 

V(q,tl£,y) =e'(«?+€2)? 2 eiq-y e -^ ^M*^. (D8) 
Finally, we invert the Fourier transform to arrive at the reduced distribution: 

P(x,f?,ff,y) =y<?qe _iqx P(q,Ci,C2 5 y) = J dqe-^- y) e^l+iDw OltoWw*. (D9) 

Convergence of this inversion is furnished by the cut-off nature of the q integration. Inserting this formula into the 
parametrization l|7UI) we arrive at a formula for the full distribution: 

7>(x,Ar u ,Ar 22 ,Ar 12 ) = J d&ti (£) ^ 2 2 AA(C 2 2 ) 6{Ar n - £ 2 l) 5(Ar 22 - £t) 

J (Jqe- iq - (x - y) eM«?+«2 2 )9 2 e -^ SWy) 9^ j 5 (Ar 12 - (T/fx )G{y)). (D10) 

I 
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